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An entanglement witness (EW) is an operator that allows 
to detect entangled states. We give necessary and sufficient 
conditions for such operators to be optimal, i.e. to detect 
entangled states in an optimal way. We show how to opti- 
mize general EW, and then we particularize our results to 
the non-decomposable ones; the latter are those that can de- 
tect positive partial transpose entangled states (PPTES). We 
also present a method to systematically construct and op- 
timize this last class of operators based on the existence of 
"edge" PPTES, i.e. states that violate the range separabil- 
ity criterion [Phys. Lett. A232, 333 (1997)] in an extreme 
manner. This method also permits the systematic construc- 
tion of non-decomposable positive maps (PM). Our results 
lead to a novel sufficient condition for entanglement in terms 
of non-decomposable EW and PM. Finally, we illustrate our 
results by constructing optimal EW acting ou H ^G"^ . 

The corresponding PM constitute the first examples of PM 
with minimal "qubit" domain, or - equivalently - minimal 
hermitian conjugate codomain. 
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I. INTRODUCTION 



Quantum entanglement which is an essence of 

many fascinating quantum mechanical effects (^-|^, is 
a very fragile phenomenon. It is usually very hard to 
create, maintain, and manipulate entangled states un- 
der laboratory conditions. In fact, any system is usu- 
ally subjected to the effects of external noise and inter- 
actions with the environment. These effects turn pure 
state entanglement into mixed state, or noisy entangle- 
ment. The separability problem, that is, the characteri- 
zation of mixed entangled states, is highly nontrivial and 
has not been accomplished so far. Even the apparently 
innocent question: Is a given state entangled and does it 
contain quantum correlations, or is it separable, and does 
not contain any quantum correlations? will, in general, 
be very hard (if not impossible!) to answer. 

Mathematically, mixed state entanglement can be de- 
scribed as follows. A density operator p > acting on a 
finite Hilbert space H — Ha 'S) Hb describing the state 
of two quantum systems A and B is called separable |^] 
(or not entangled) if it can be written as a convex com- 
bination of product vectors; that is, in the form 



P = ^Pk\ek,fk){ek,fk\, 



(1) 



where pk > 0, and \ek, fk) = \ek)A ® |/fe)s are product 
vectors. Conversely, p is nonseparable (or entangled) if 
it cannot be written in this form. Physically, a state de- 
scribed by a separable (nonseparable) density operator 
p can always (never) be prepared locally. Most of the 
applications in quantum information are based on the 
nonlocal properties of quantum mechanics, and 
therefore on nonseparable states. Thus, a criterion to 
determine whether a given density operator is nonsepa- 
rable, i.e. useful for quantum information purposes, or 
not is of crucial importance. On the other hand, PPTES 
are objects of special interest since they represent so- 
called bound entangled states, and therefore provide an 
evidence of irreversibility in quantum information pro- 
cessing Q. 

For low dimensional systems there exist oper- 

ationally simple necessary and sufficient conditions for 
separability In fact, in = (T^ 0(7^ and H ^(F'^ (^(F^ 
the Peres-Horodecki criterion [ p"oHTl] ] establishes that p 
is separable iff its partial transpose is positive. Partial 
transpose means a transpose with respect to one of the 
subsystems ||l2[. For higher dimensional systems all op- 
erators with non-positive partial transposition are en- 
tangled. However, there exist positive partial transpose 
entangled states (PPTES) ||l|,|l|l. Thus, the separabil- 
ity problem reduces to finding whether density opera- 
tors with positive partial transpose are separable or not 

In the recent years there has been a growing effort in 
searching for necessary and sufficient separability criteria 
and checks which would be operationally simple |]l5| , |l6| . 
Several necessary [0jl^] or sufficient ]l3| , p^| -pll conditions 
for separability are known. A particularly interesting 
necessary condition is given by the so-called range cri- 
terion |l^. According to this criterion, if the state p 
acting on a finite dimensional Hilbert space is separable 
then there must exist a set of product vectors {\ek, fk)} 
that spans the range R{p) such that the set of partial 
complex conjugated product states {\ek,fk)} spans the 
range of the partial transpose of p with respect to the 
second system, i.e., p^^ . Among the PPTES that vi- 
olate this criterion there are particular states with the 
property that if one subtracts a projector onto a prod- 
uct vector from them, the resulting operator is no longer 
a PPTES (20|j2^. In this sense, these states fie in the 
edge between PPTES and entangled states with non- 
positive partial transposition, and therefore we will call 
them "edge" PPTES. The analysis of the range of den- 
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sity operators initiated in Ref. has turned out to be 
very fruitful. In particular, it has led to an algorithm for 
the optimal decomposition of mixed states into a separa- 
ble and an inseparable part [^-^, and to a systematic 
method of constructing examples of PPTES using unex- 
tendible product bases . For low rank operators it 

has allowed to show that one can reduce the separability 
problem to the one of determining the roots of certain 
complex polynomial equations . 

From a different point of view, a very general approach 
to analyze the separability problem is based on the so- 
called entanglement witnesses (EW) and positive maps 
(PM) Entanglement witnesses are operators 

that detect the presence of entanglement. Starting from 
these operators one can define PM's that also detect 
entanglement. An example of a PM is precisely partial 
transposition |]l0| , ^ , ^ . The importance of EW stems 
from the fact that a given operator is separable iff there 
exists an EW that detects it Thus, if one was able 
to construct all possible EW (or PM) one would have 
solve the problem of separability. Unfortunately, it is 
not known how to construct EW that detect PPTES in 
general. The only result in this direction so far has been 
given in Ref. |25|, although some preliminary results ex- 
ist in the mathematical literature Starting from a 
PPTES fulfilling certain properties (related to the exis- 
tence of uncxtendible basis of product vectors ||l^]), it 
has been shown how to construct an EW (and the corre- 
sponding PM) that detects it. Perhaps, one of the most 
interesting goals regarding the separability problem is to 
develop a constructive and operational approach using 
EW and PM that allows us to detect mixed entangle- 
ment. 

In this paper we realize this goal partially: we intro- 
duce a powerful technique to construct EW and PM that, 
among other things, allows us to study the separability of 
certain density operators . In particular, we show how to 
construct optimal EW; that is, operators that detect the 
presence of entanglement in an optimal way. We specif- 
ically concentrate on non-decomposable EW, which are 
those that detect the presence of PPTES. Furthermore, 
we present a way of constructing optimal EW for edge 
PPTES. Our method generalizes the one introduced by 
Terhal |2^] to the case in which there are no uncxtendible 
basis of product vectors. When combined with our pre- 
vious results 1^,^ regarding subtracting product vec- 
tors from PPTES, the construction of non-decomposable 
optimal EW starting from "edge" PPTES gives rise to 
a novel sufficient criterion for non-separability of gen- 
eral density operators with positive partial transposi- 
tion. We illustrate our method by constructing optimal 
EW that detect some known examples of PPTES 
in H = ^C^. The corresponding PM constitute the 
first examples of PM with minimal "qubit" domain, or - 
equivalently - minimal hermitian conjugate codomain. 

This paper is organized as follows. In Section II we 
review the definition of EW and fix some notation. In 
Section III we study general EW. We define optimal wit- 



nesses and find a criterion to decide whether an EW is 
optimal or not. In Section IV we restrict the results of 
Section III to non-decomposable EW. In particular, we 
show how to optimize them by subtracting decomposable 
operators. In Section V we give an explicit method to 
optimize both, general and non-decomposable EW. We 
also show how to construct non-decomposable EW, and 
that this leads to a sufficient criterion of non-separability. 
The construction and optimization is based on the use of 
"edge" PPTES. In Section VI we extend our results to 
positive maps. In Section VII we illustrate our methods 
and results starting from the examples of PPTES given 
in Ref. |jl^. The paper also contains two appendices. 
In Appendix A we describe in detail a method to check 
whether an EW is optimal or not. In Appendix B we 
discuss separately some important properties of the edge 
PPTES, and show that they provide a canonical decom- 
position of mixed states with positive partial transpose. 



II. DEFINITIONS AND NOTATION 



We say that an operator W 
Ha ® Hb is an EW if f^^: 



acting on H 



(I) (e, /|VF|e, /) > for aU product vectors |e, /); 

(II) has at least one negative eigenvalue (i.e. is not pos- 

itive) ; 

(III) tr(VK) = 1. 

The first property (I) imphes that {p)w = tr{Wp) > 
for all p separable. Thus, if we have {p)w < for some 
p > 0, then p is nonseparable. In that case we say that 
W detects p. The second one (II) implies that every 
EW detects something, since in particular it detects the 
projector on the subspace corresponding to the negative 
eigenvalues of W. The third property (III) is just nor- 
malization condition that we need in order to compare 
the action of different EW |§. 

In this paper we will denote by K{p) and R{p) the ker- 
nel and range of p, respectively. The partial transposition 
of an operator X will be denoted by X'^ p^ , |3l[ |. On the 
other hand, we will encounter several kinds of operators 
(EW, positive operators, decomposable operators, etc) 
and vectors. In order to help to identify the kind of op- 
erators and vectors we use, and not to overwhelm the 
reader by specifying at each point their properties, we 
will use the following notation: 

• W will denote an EW. 

• P,Q will denote positive operators. Unless speci- 
fied they will have unit trace [tr(P) = tr((5) = 1]. 

• D will denote a decomposable operator. That is, 
D = aP + bQ^ , where a, 6 > 0. Unless stated, all 
decomposable operators that we use will have unit 
trace (i.e., 6 = 1 — a). 
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• p will denote a positive operator (not necessarily of 
trace 1). 

• |e, /) will denote product vectors with |e) G Ha 
and I/) G Hb- Unless especified, they will be nor- 
malized. 



III. GENERAL ENTANGLEMENT WITNESSES 

In this Section we first give some definitions directly 
related to EW. Then we introduce the concept of opti- 
mal EW. We derive a criterion to determine when an 
EW is optimal. This criterion will serve us to find an 
optimization procedure for these operators. 

A. Definitions 

Given an EW, W, we define: 

• Dw = {p > 0, such that {p)w < 0}; that is, the 
set of operators detected by W. 

• Finer: Given two EW, Wi and W2, we say that W2 
is finer than Wi , if Q ; that is, if all the 
operators detected by Wi are also detected by W2. 

• Optimal entanglement witness (OEW): We say 
that W is an OEW if there exist no other EW which 
is finer. 

• Pw = {|e,/) e H, such that {eJ\W\eJ) = 0}; 
that is, the set of product vectors on which W van- 
ishes. As we will show, these vectors are closely 
related to the optimality property. 

Note the important role that the vectors in Pw play 
regarding entanglement (for a method to determine Pw 
in practice, see Appendix A). If we have an EW, W , 
which detects a given operator p, then the operator p' = 
p + Pw where 



Pw = y^Pfc|efc,/fc)(efc,/fc| 



(2) 



with pk > 0, and \ek,fk) G Pw is also detected by W. In 
fact, this means that any operator of the form ^ is in 
the border between separable states and non-separable 
states, in the sense that if we add an arbitrarily small 
amount of p to it we obtain a non-separable state. Thus, 
the structure of the sets Pw characterizes the border be- 
tween separable and non-separable states. In fact, from 
the results of this Section it will become clear that we 
can restrict ourselves to the structure of the sets of Pw 
corresponding to OEW's. 



B. Optimal entanglement witnesses 

According to Ref. p is nonseparable iff there exists 
an EW which detects it. Obviously, we can restrict our- 
selves to the study of OEW. For that, we need criteria to 
determine when an EW is optimal. In this subsection we 
will derive a necessary and sufficient condition for this to 
happen (Theorem 1 below). In order to do that, we first 
have to introduce some results that tell us under which 
conditions an EW is finer than another one. 

Lemma 1: Let W2 be finer than Wi and 



A; 



inf 



{P 



IIW2 



(P 



(3) 



Then we have: 



(i) If {p)w, = then {p)w, < 0. 

(ii) If {p)wi < 0, then {p)w2 < {p)wi- 

(iii) If {p)wi > then X{p)wi > {p)w2- 

(iv) A > 1. In particular, A = 1 iff I^i = W2. 

Proof: Since W2 is finer than Wi we will use the fact 
that for all p > such that {p)wi < then {p)w2 < 0. 

(i) Let us assume that {p)w2 > 0- Then we take any 
Pi € T^Wi so that for all a; > 0, < p{x) = pi -\- xp £ 
T^Wi- But for sufficiently large x we have that {p{x))w2 
is positive, which cannot be since then p(x) ^'Dw2 ■ 

(ii) We define p = p + \ {p)wi\^ > 0. We have that 
{p)wi = 0. Using (i) we have that > {p)w2 + l(p)wi |- 

(iii) We take pi £ and define p — {p)wiPi + 
|(Pi)m^iIp > 0, so that {p)wi = 0. Using (i) we have 
\{pi)wi\{p)w2 < \{pi}w2\{p)wi- Dividing both sides by 

> (p)wi > we obtain 



(P> 



W2 



(P) 



< 



Wi 



{P 



1IW2 



{P 



l)Wi 



(4) 



Taking the infimum with respect to pi G Dwi the rhs 
of this equation we obtain the desired result. 

(iv) From (ii) immediately follows that A > 1. On the 
other hand, we just have to prove that if A = 1 then Wi = 
W2 (the only if part is trivial). If A = 1, using (i) and 
(in) we have that {pv)wi > {Pv)w2 for all p„ = |e, /)(e, /| 
projector on a product vector. Since tr(V7i) = tr(W2) we 
must have tr[(M^i — W2)pv] = for all p„, since we can 
always find a product basis in which we can take the 
trace. But now, for any given p > we can define p{x) — 
p+xl such that for large enough x, p{x) is separable p8[ . 
In that case we have {p{x))wi = {p{x))w2 which implies 
that {p)w^ = (p)w2, i-e- VFi = Wa- ^ 

Corollary 1: Dw^ = Dw2 iff Wi = W2. 

Proof: We just have to prove the only if part. For that, 
we define A as in (0). On the other hand, defining 



A = inf 

P2&DW2 



(P 



2)Wi 



(P2) 



14^2 



(5) 
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we have that A > 1 since Wi is finer than W2 (Lemma 
l(iv)). Equivalently, 



1 > sup 



{P 



1)W2 



{P 



1/Wi 



> A > 1, 



(6) 



where for the last inequaUty we have used that W2 is 
finer than Wi. Now, since A = 1 we have that Wi = W2 
according to Lemma l(iv). □ 

Next, we introduce one of the basic results of this pa- 
per. It basically tell us that EW is finer than another 
one if they differ by a positive operator. That is, if we 
have an EW and we want to find another one which is 
finer, we have to subtract a positive operator. 

Lemma 2: W2 is finer than Wi iff there exists a P 
and 1 > e > such that Wi = {1 - e)W2 + eP. 

Proof: (If) For all p G Dw^ we have that > {p)wi — 
(l—e)(p)vF2+e(p)p which implies {p)w2 < and therefore 
p e Dw2- (Only if) We define A as in (||). Using Lemma 
l(iv) we have A > 1. First, if A = 1 then according to 
Lemma l(iv) we have Wi = W2 (i.e., e = 0). For A > 1, 
we define P= (A- l)-i(AI¥i - W2) and e = 1-1/A> 0. 
We have that Wi — {l—e)W2+f-P, so that it only remains 
to be shown that P > 0. But this follows from Lemma 



A = inf 



PieDv 



(pi> 



l(i-iii) and the definition of A, 

□ 

The previous lemma provides us with a way of deter- 
mining when an EW is finer than another one. With this 
result, we are now at the position of fully characterizing 
OEW. 

Theorem 1: 14^ is optimal iff for all P and e > 0, 
W = {l + e)W - eP is not an EW [does not fulfiU (I)]. 

Proof: (If) According to Lemma 2, there is no EW 
which is finer than W, and therefore W is optimal. (Only 
if) If W is an EW, then according to Lemma 2 is not 
optimal. □ 

The previous theorem tells us that W is optimal iff 
when we subtract any positive operator from it, the re- 
sulting operator is not positive on product vectors. This 
result is not very practical because of two reasons: (1) 
for a given P it is typically very hard to check whether 
there exists some e > such that W — eP is positive 
on all product vectors; (2) it may be difficult to find a 
particular P that can be subtracted from W among all 
possible positive operators. In Appendix A we show how 
to circumvent these two drawbacks in practice: we give 
a simple criterion to determine when a given P can be 
subtracted from W. This allows us to determine which 
are the positive operators which can be subtracted from 
a given EW. 

In the rest of this subsection we will present some sim- 
ple results related to these two questions. First, it is clear 
that not every positive operator P can be subtracted 
from an EW, W. In particular, the following lemma tells 
us that it must vanish on Pw- 

Lemma 3: If PPw 7^ then P cannot be subtracted 
from W. 



Proof: There exists some |eo,/o) G Pw such that 
(eoj /o|-P|eoj /o) > 0. Substituting this product vector 
in the condition I for any — eP we see that the in- 
equality is not fulfilled for any e > 0, i.e. P cannot be 
subtracted. □ 

Corollary 2: If Pw spans H then W is optimal. 

Note that, as announced at the beginning of this Sec- 
tion, the set Pw plays an important role in determining 
the properties of the separable states which lie on the 
border with the entangled states. We see here, that this 
set also plays an important role in determining whether 
an EW is optimal or not. 

On the other hand, in order to check whether a given 
operator P can be subtracted or not from W, one has to 
check whether there exists some e > such that (e, /| VF— 
eP|e,/) > for all |e,/). The following lemma gives an 
alternative way to do this. In fact, it gives a necessary 
and sufficient criterion for an EW to be optimal. For a 
given |e) £ Ha, we will denote by We = {e\W\e). 

Lemma 4: W is optimal iff for all |^) orthogonal to 
Pw 



inf r(*i 



= 0. 



(7) 



Proof: (If) Let us assume that W is not optimal; that 
is, there exists W ^ W, finer than W. Then, accord- 
ing to Lemma 2 we have that there exists eg > and 
P > such that W' = {W - eoP)/(l - eo). Imposing 
that W is positive on product vectors (i.e. > 

for all |e) € Ha) we obtain < {e\W - eQP\e) < 
We — eoA3i(e|^)(4'|e), where l^*) is any eigenstate of 
P with nonzero eigenvalue Aij,. According to Ref. [ pO| , 
this last operator is positive iff both: (i) (e|^) is in the 
range of (e|iy|e), which imposes that |^) is orthogonal 

to Pw; (ii) Avpeo < [(*|e)I^-i(e|*)]" , which imposes 
that e > A^eo > for that given f^f). (Only if) Let 
us assume that there exists some j^") orthogonal to Pw 
such that e > 0. Then, using the same arguments one 
can show that W' = {W ~ e|*)(*|)/(l - e) ^ W is an 
EW. According to Lemma 2, W' is finer than W, so that 
W is not optimal. □. 



C. Decomposable entanglement witnesses 

There exists a class of EW which is very simple to char- 
acterize, namely the decomposable entanglement wit- 
nesses (d-EW) Us). Those are EW that can be written 
in the form 



W ^aP+{l- a)Q^, 



(8) 



where a G [0, 1]. As it is well known (see next section), 
these EW cannot detect PPTES. In any case, for the sake 
of completeness, we will give some simple properties of 
optimal d-EW. 
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Theorem 2: Given a d-EW, W, if it is optimal then 
it can be written a.s W = Q'^ , where Q > contains no 
product vector in its range. 

Proof: Since W is decomposable, it can be written 
as W =^ aP + {1 - a)Q^. W (x W - aP is also a 
witness, which according to Lemma 2 is finer than W, 
and therefore W is not optimal. On the other hand, 
if |e, /) € R{Q) then for some A > we have that 
W (x {Q — A|e, /) (e, /l)"^ is finer than W, and therefore 
this last is not optimal. □ 

This previous result can be slightly generalized as fol- 
lows: 

Theorem 2': Given a d-EW, W, if it is optimal then 
it can be written as 14^ = , where Q > and there is 
no operator P S R{Q) such that P^ > 0. 

Proof: Is the same as in previous theorem. □ 

Corollary 3: Given a d-EW, W, if it is optimal then 
W"^ is not an EW [does not fulfill (II)]. 

Proof: Using Theorem 2 we have that W = with 
g > 0. Then W'^ = Q > 0, which docs not satisfy 
property (ii). □ 



IV. NON-DECOMPOSABLE ENTANGLEMENT 
WITNESSES 

In the previous section we have been concerned with 
EW in general. As mentioned above, when studying sep- 
arability we just have to consider those EW that can 
detect PPTES. In order to characterize them, one de- 
fines non-decomposable witnesses (nd-EW) as those EW 
which cannot be written in the form (^) |2|]. This Sec- 
tion is devoted to this kind of witnesses. The importance 
of nd-EW in order to detect PPTES is reflected in the 
following 

Theorem 3: An EW is non-decomposable iff it de- 
tects PPTES. 

Proof: (If) Let us assume that the EW is decompos- 
able. Then it cannot detect PPT, since if p, > we 
ha.veti[{aP+{l-a)Q'^)p] = atT{Pp) + {l-a)tr{Qp^) > 
0. (Only if) The set of decomposable witnesses is convex 
and closed, and VF, as a set containing one point, is a 
closed convex set itself. Thus, from Hahn-Banach the- 
orem it follows that there exists an operator p such 
that: (i) tr[p(aP + (1 - a)Q^)] > for all P,Q > 0, 
a G [0, 1]; (ii) tr{pW) < 0. From (i), taking a = 1 we in- 
fer that p > 0; on the other hand, taking a = we obtain 
that ti[p'^Q] > for all Q > 0, and therefore p'^ > 0. 
Thus, W detects p which is a PPTES. □ 

Corollary 4: Given an operator D, it is decomposable 



iff tr(L>p) > for aU p, p^ > 0. 



A. Definitions 

In this Subsection we introduce some definitions which 
are parallel to those given in the previous Section. Given 



a nd-EW, W, we define: 

• dw = {p > 0, such that p^ > and {p)w < 0}; 
that is, the set of PPT operators detected by W. 

• Non-decomposable-finer (nd-finer): Given two nd- 
EW, Wi and W2, we say that W2 is nd-finer than 
Wi, if dvFi ^ dw2'i that is, if all the operators de- 
tected by Wi are also detected by W2- 

• Non-decomposable optimal entanglement witness 
(nd-OEW): We say that W is an nd-OEW if there 
exist no other nd-EW which is nd-finer. 

'PW = {|e,/) e H, such that {eJ\W\eJ) = 0}; 
that is, the product vectors on which W vanishes. 

Note again the important role that the vectors in pw 
play regarding PPTES. If we have a nd-EW, W, which 
detects a given PPTES p, then the operator p' ~ p + Pw 
where pw has the form (||) with pk > 0, and \ek,fk) G 
pw also describes a PPTES. Thus, any operator of the 
form (^ lies in the border between separable states and 
PPTES. 



B. Optimal non— decomposable entanglement witness 

The goal of this section is to find a necessary and suf- 
ficient condition for a nd-EW to be optimal. We start 
by proving a similar result to the one given in Lemma 1, 
but for nd-EW: 

Lemma lb: Let W2 be nd-finer than Wi, 



A= inf 



tr(W^2Pi) 



tr(VFiPi) 



(9) 



and now both, p, p^ > 0. Then we have have (i-iv) as in 
Lemma 1. 

Proof: The proof is basically the same as in Lemma 1 
and will be omitted here. 

Corollary lb: Given two nd-EW, M^i,2, then dwi — 
dw2^SWi=W2. 

Proof: The proof is basically the same as Corollary 1 
and will be omitted here. 

Lemma 2b: Given two nd-EW, Wi,2, W2 is nd-finer 
than Wi iff there exists a decomposable operator D and 
1 > e > such that Wi ^ {1 ~ e)W2 + eD. 

Proof: (If) Given any p, p^ > 0, we have that ii p G 

dwi then > {p)wi = (1 - e)(p)w2 + > {p)w2, 

where in the last inequality we have used that {p)d > 
since D is decomposable (see Corollary 4). Therefore 
p € dw2- (Only if) We define A as in (||), so that A > 1 
according to Lemma Ib(iv). If A = 1 we have Wi = W2- 
If A > 1 we define D = (A - l)-^{\Wi - W2) and e 
1 - 1/A. We have that VFi = (1 - 6)1^2 + £-0, so that it 
only remains to be shown that D is decomposable. But 
from Lemma Ib(i-iii) and the definition of A it follows 
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that {p) D ^ for all p, > 0. Using Corollary 4 we 
then have that D is decomposable. □. 

Now we are able to fully characterize nd-OEW. 

Theorem lb: Given an nd-EW, W, it is nd-optimal 
iff for all decomposable operators D and e > 0, W = 
(1 + e)W^ - eL* is not an EW [does not fulfill (I)]. 

Proof: Is the same as for Theorem 1. □ 

Theorems 1 and lb allow us to relate OEW and nd- 
OEW. In this way we can directly translate the results 
for general OEW to nd-OEW. We have 

Theorem 4: Given a nd-EW, W, W is a nd-OEW 
iff both W and are OEW. 

Proof: (If) Let us assume that W is not a nd-OEW. 
Then, according to Theorem lb there exists e > and 
a decomposable operator D such that W = {1 + e)W — 
eD is a nd-EW. We can write D = aP + {1 - a)Q'^ , 
with a £ [0, 1]. If a 7^ 0, then T^i = (1 + ae)W - aeP 
fulfills {eJ\Wi\e,f) > {eJ\W'\e,f) > 0, and therefore, 
according to Lemma 2, W is not optimal. If a 7^ 1 then 
W2 = [H-(l-a)e]VF^-(l-a)eQ fulfills (e, /| W^ale, /) > 
{eJ\{W')^\eJ) > 0, i.e. is an EW and therefore 
is not optimal. (Only if) According to Theorem lb, if 
W is nd-optimal then for all D = aP + (1 — a)Q^ , with 
a € [0, 1], and all e > we have that W = {1- e)W-eD 
does not satisfy (I). Taking a = 1 we have for all P and 
e > 0, H^i = (1 - e)W - eP does not fulfill (I), and 
therefore (Theorem 1) W is optimal; analogously, taking 
a = we have that is optimal also. □ 

Corollary 5: W is a nd-OEW iff is an nd-OEW. 



V. OPTIMIZATION 

In this Section we give a procedure to optimize EW 
which is based on the results of the previous Sections. 



A. Optimization of general entanglement witnesses 

Our method is based in the following lemma. It tells 
us how much we can subtract from an EW. Here we will 
denote by We — {e\W\e) and Pe = (e|P|e) where |e) G 
Ha, by [. . .]min the minimum eigenvalue, and by [. . .]max 
the maximum eigenvalue. On the other hand, X^^/^ will 
denote the square root of the pseudoinverse of X Q . 

Lemma 5: If there exists some P such that PPw — 
and 



Aq = inf 

\e)eHA 



p-^/^WeP^^^^ 



(10) 



Proof: Let us find out for which values of A > 0, W'{\) 
is an EW. We have to impose condition (I), which can 
be written as (e|W(A)|e) > 0, i.e. 



W, - AR > 0. 



— 1/2 

Multiplying by Pe on the right and left of this equa- 



(12) 



tion we obtain Pe ^^'^WePe '''^ > A, which immediately 
gives that A < Aq given in the first part of Eq. (|l^). 

— 1/2 

On the other hand, multiplying by We on the right 
and left of Eq. (|l|) we obtain W^^^^ PeW;r'^^'^ < 1/A, 
which immediately gives that A < Aq given in the second 
equality of Eq. (0). □ 

Lemma 5 provides us with a direct method to optimize 
EW by subtracting positive operators for which the ele- 
ments of Pw are contained in their kernels. The method 
thus consists of: (1) determining Pw', (2) choosing an op- 
erator P so that PPw — and determining A using (p^); 
(3) if A 7^ then we subtract the operator P according to 
Lemma 5. Continuing in the same vein we will reach an 
OEW. In Appendix A we show how to accomplish steps 
(1) and (2) in practice. 



1/2 



B. Optimization of non— decomposable entanglement 
witnesses 

For nd-EW we have the following generalization of 
Lemma 5: 

Lemma 5b: Given a nd-EW, W, if there exists some 
decomposable operator D such that Dpw = and 



Aq = inf 



D-^/^WeD-^/^ 



mm 
-1 



(13) 



sup 



We^^'^DeW-^l'^ 



> 0. 



then 



l^'(A) = (M^ - A7:>)/(1 - A) 



(14) 



with A > is a nd-EW iff A < Aq. 

Proof: Is the same as for Lemma 5. 

With the help of Lemma 5b we can optimize nd-EW by 
subtracting decomposable operators as follows: (1) deter- 
mining Pw and pwT ; (2) choosing P, Q so that Ppw = 
and QpwT = 0, building D = aP + (1 — a)Q^ with 
a £ [0, 1], and determining Aq using (|l^); (3) if Aq 7^ 
then we subtract the operator D according to Lemma 5b. 



sup 



W^^/^PeWe^^^ 



> 0. 



C. Detectors of "edge" PPTES 



then 

W'{\) = {W - XP)/{1- X) 
with A > is an EW iff A < Aq. 



In the previous subsections we have have given two op- 
^21) timization procedures. In both of them, starting from a 
general EW one can obtain one which is optimal (or nd- 
optimal) . It may well happen that the EW found in this 
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way is non-decomposable even though the original one 
was decomposable. To check that one simply has to use 
Corollary 3; that is, check whether is an EW or not. 
In case it is, then the OEW W is non-decomposable. 
However, nothing guarantees that the final EW is non- 
decomposable if the original one is not. In this subsection 
we describe a general method to construct nd-EW us- 
ing the optimization procedures introduced earlier. This 
method generalizes the one presented in Ref . |25j . 

We are going to use the results presented in Ref. 
There, we have already used and discussed the 
"edge" PPTES, without naming them, however. Let us 
now introduce the following definition: 

Definition: [see Ref. [Q] A PPTES S is an "edge" 
PPTES if for all product vector |e, /) and e > 0, (5 - 
e|e,/)(e,/| is not a PPTES. 

This definition implies that that the "edge" states lie 
on the boundary between the PPTES and entangled 
states with non-positive partial transpose. In this sub- 
section we will show how, out of an "edge" PPTES, we 
can construct a nd-OEW that detects it. As we men- 
tioned in the introduction, "edge" PPTES are of special 
importance. In particular, they allow to provide a canon- 
ical form to write PPTES in arbitrary Hilbert spaces. 
For these reasons, some of the properties of the "edge" 
PPTES are discussed in Appendix B. 

In order to check whether a PPTES S is an "edge" 
PPTES we can use the range criterion ||l3| (see also ) . 
That is, S is an "edge" PPTES iff for aU |e, /) G R{6), 

Let 6 be an "edge" PPTES, and let us denote by Pi 
the projector onto K{S) and by Qi the projector onto 
KiS'^). We define 

Ws - a(Fi + Qj), (15) 

where a = l/tr(Pi -I- Qi). Let us also define 

ei= inf (e,/|M^,|e,/). (16) 

Then we have 

Lemma 6: Given an "edge" PPTES S, then Wi cx 
Ws — eil is a nd-EW, where ci and Ws are defined in 
(^,^, respectively. 

Proof: We have that (e, f\Ws\e, /) = a((e, f\Pi\e, /) + 
(e,/*|Qi|e,/*)) > 0. This quantity is zero iff 
(e, /|Pi|e,/) = (e, /*|(5i|e, /*) =0. But this is not possi- 
ble since S is an "edge" PPTES. Thus, (e, f\Ws\e, /) > 
for all |e,/). Defining ci as in (^6|), and taking into ac- 
count that (e, /iPV^je, /) is a continuous function of (the 
coefficients of) |e,/) and that the set in which we are 
taken the infimum is compact, we obtain ei > 0. Then we 
obviously have that Wi fulfills properties (I) and (III). On 
the other hand, (15)1^1 oc a{{6)p^ +(^"^)Qi)~£i < 0: since 
PiS = QiS'^ = 0. Thus, Wi detects a PPTES, and there- 
fore, according to Theorem 3 is non-decomposable. □ 

Note that Lemma 6 provides an important generaliza- 
tion of the method of Terhal [E5|, based on the use of 



unextendible product bases jlj]. Our method works in 
Hilbert spaces of arbitrary dimensions, and in particu- 
lar when dim(i?yi) = 2 (in 2 x dimensional systems) 
for which unextendible product basis do not exist. By 
combining Lemma 6 and the optimization procedure in- 
troduced earlier, we obtain a way of creating nd-OEW. 
Once we have Wi we find pwi and py^^T . We denote by 
P2 and Q2 the projector operators orthogonal to these 
two sets, respectively. 



£2 



inf 



{eJ\Wi\eJ) 



J\P2 



(17) 



and W2 oc Wi — £2(^2 + Qi")- According to Lemma 2b 
we have that W2 is nd-finer than Wi . Now we can define 
PW2 1 Pw^ T P^tQ^ and W3 in the same way, and continue 
in this vein until for some fc, = 0. If Wk is not yet 
optimal, we still have to find other projectors such that 
we can optimize as explained in the previous subsections. 

In Section VII we illustrate this method with a family 
of edge PPTES from Ref. [|3|. In fact, as we will mention 
in that Section, we have checked that the optimization 
method typically works as well by starting with three 
random vectors, and following a similar procedure to the 
one indicated here. This means that in our construction 
method we do not need in practice to start from a "edge" 
PPTES. 



D. Sufficient condition for PPTES 

In this subsection we use the results derived in the 
previous one to construct a sufficient criterion for non- 
separability of PPTES. As shown in Ref. [|o|j2l), given an 
operator p > 0, with p"^ > 0, we can always decompose 
it in the form 

P^Ps + S, (18) 

where ps is separable and S is an "edge" PPTES. More 
details concerning this decomposition, and in particular 
its canonical optimal form are presented in Appendix B. 
In this section we use this decomposition together with 
the following 

Lemma 7: Given a non-separable operator p = p^ + S, 
where > is separable then for all EW, W, such that 
{p)w < we have that {d)w < 0. 

Proof: Obvious from the definition of EW. □ 
Lemma 7 tells us that if p is non-separable, then there 
must exist some EW that detect both S and p. Actually, 
it is clear that there must exist an OEW with that prop- 
erty. In particular, if p^ > 0, it must be a nd-OEW. 
In the previous subsection we have shown how to build 
them out of "edge" PPTES. Thus, given p we can always 
decompose it in the form (p^), construct an OEW that 
detects 5 and check whether it detects p. In that case, 
we will have that p is non-separable. Thus, this provides 
a sufficient criterion for non-separability. 
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We stress the fact that for PPTES only a special class 
of states, namely the class of "edge" PPTES, is responsi- 
ble for the entanglement properties. In fact, one should 
stress that very many of the examples of PPTES dis- 
cussed so far in the literature bel ong to the class of "edge" 
PPTES: the 2 (g) 4 family from the n » n states ob- 
tained via unextendible product basis construction p4[ , 
the 3 ® 3 states obtained via the chess-board method 
|Q(b), and projections of continuous variable PPTES 
onto finite dimensional subspaces p3] (c) . 



VI. POSITIVE MAPS 

It is known that PM allow for necessary and suffi- 
cient conditions for separability (or, equivalently, entan- 
glement) of bipartite mixed states PM's have been 
also applied in the context of distillation of entangle- 
ment 1 35 1 and information theoretic analysis of separa- 
bility 1 3f ] . In this Section we will use the isomorphism 
between operators and linear maps to extend the proper- 
ties derived for witnesses to PM . We will first review 
some of the definitions and properties of linear maps. 

Let us consider a linear map £ : B{Ha) B{Hc)- 
We say that £ is positive if for all Y € B{Ha) posi- 
tive, £iY) > 0. One can extend a linear map as fol- 
lows. Given £ : B{Ha) B{Hc), we define its extension 
£®Ib: B{Ha) ® B{Hb) ^ B{Hc) ® B{Hb) according 
to £(^lc{Y.^Y^®Zi) ^Y.^^{yi)®Z,,^\veYeY, e B{Ha) 
and Zi e B{Hb)- A linear map is completely positive if 
all extensions are positive. The classification and charac- 
terization of positive (but not completely positive) maps 
is an open question (see, e.g. Ref. p^ , |29| ]). 

An example of positive (but not completely positive) 
map is transposition (in a given basis Oa)] that is, the 
map £t such that £t{X) — ■ The corresponding ex- 
tension is the partial transposition . A map £ is called 
decomposable if it can be written a,s £ ^ £i + £2 ■ £t, 
where £"1^2 are completely positive. 

One can relate linear maps with linear operators in 
the following way. We will assume cIa = dim(_ffyi) < 
dim(i?c'), but one can otherwise exchange Ha by He in 
what follows. Given X € B{Ha ® He) and an orthonor- 
mal basis Oa = {\k)}iti in Ha, we define the linear map 
£x ■ B{Ha) — > B{Hc) according to 



£(y) -tr^(x^^y). 



(19) 



for all Y e B{Ha), where tr^i denotes the trace in Ha 
and the partial transpose is taken in the basis Oa- Sim- 
ilarly, given a linear map we can always find an operator 
X such that (noh is fulfilled. For instance, if we choose 
T = (!«')(* I )*7^where 



1*) 



(E> \k)c 



(20) 



then the corresponding map £t is precisely the transpo- 
sition in the basis OA- 
Given a linear map £x , one can easily show the follow- 
ing relations: (a) £x is completely positive iff X > 0; (b) 
£x is positive but not completely positive iff X is an EW 
[except for the normalization condition (III)]; (c) £x is 
decomposable iff X is decomposable. Thus, the problem 
of studying and classifying PM is very much related to 
the one of EW. Furthermore, PM can be also used to 
detect entanglement [|ll[. Let us consider the extension 
£x = £x®l ■- B{Ha) ® B{Hb) ^ B{Hc) ® B{Hb), 
where we take (Ib = dim(_ffB) = dim(iJc'). Then we 
have that given p € B{Ha <8) Hb), 



{p)x = {^\£x{pm, 



where 



\^)=Y.\k)c®\k)B- 



(21) 



(22) 



k=l 



Thus, if an EW, W, detects p, then £w{p) is a non- 
positive operator. Consequently, p > is entangled iff 
there exists a PM such that acting on p gives a non- 
positive operator. In that case we say that the PM "de- 
tects" p- Actually, PM are "more efficient" in detecting 
entanglement than EW. The reason is that it may hap- 
pen that £x{p) is non-positive but still {p)x > 0. 

It is convenient to define finer and optimal PM as for 
EW. That is, given two PM, £1^2, we say that £2 is finer 
than £1 if it detects more. We say that a PM, £, is 
optimal if there exists no one that is finer. In the same 
way we can define nd-finer and nd-optimal. 

The results presented in the previous sections can be 
directly translated to PM given the following fact. 

Lemma 8: If W2 is finer (nd-finer) than Wi then £w2 
is finer (nd-finer) than £wi - 

Proof: Using Lemma 2 we can write Wi = (1 — e) W2 + 
eP. According to (|l^) we have that £wi = (1 — <^)£w2 + 
€£p. Since £p{p) > for all p > 0, we have that £w2 is 
finer than £wi - Using Lemma 2b we can also prove that 
it is nd-finer. □ 

^From this lemma it follows that optimizing EW im- 
plies optimizing PM. In fact, the constructions that we 
have given in the Section VG can be viewed as ways 
of constructing non-decomposable PM. In fact, since 
the method works for dim(iJ^) = 2, the resulting PM 
£ : B{Ha) — > B{Hc) has a minimal "qubit" domain, or 
- equivalently - minimal hermitian conjugate codomain. 
Up to our knowledge, our method is the first one that 
permits to construct non-decomposable PM with these 
characteristics. 



VII. ILLUSTRATION 

In this section we explicitly give construct a nd-OEW 
out of edge PPTES. We use, as a starting point, the 
family of PPTES introduced in |l|]). 
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A. Family of "edge" PPTES 

We consider Ha = and Hb = <I'^, and denote 
by {|fc)}fc=o — ^i-S) an orthonormal basis in these 
spaces, respectively. Most of the time we will write the 
operators in those bases; that is, as matrices. For oper- 
ators acting in Ha ® Hb we will always use the follow- 
ing order {|0,0), |0, 1), . . . , |1, 0), |1, 1), . . .}. On the other 
hand, all partial transposes will be taken with respect to 
Hb- 

We consider the following family of positive operators 
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where 6 E [0, 1]. For 6 = 0, 1 those states are separable, 
whereas for < 6 < 1, is an "edge" PPTES. This can 
be shown by checking directly that they violate the range 
criterion of Ref. i.e. the definition given in Section 
IVC. 

If we take the partial transpose in the basis 
the density operators ph have the property that — 

UBPbUg with Ub = {crx)o3 © {(Tx)i2- Here, the subscript 
ij denotes the subspace, Hsij C Hb spanned by {\i), \j)} 
and ax is one of the Pauli-operators. Note that Ub is a 
real unitary operator acting only on Hb- This immedi- 
ately implies that 

pI = Pb, (24) 

where pb = VBPbVl and Vb = 1^2^ + iax)m ® (1 + 
i<7x)i2]- We will use the property (gj) to simplify the 
problem of constructing the nd-OEW. Thus, we will con- 
centrate from now on the operators pb WA - Obviously, 
Pb is an "edge" PPTES for 1 > 6 > 0. 

The projector onto the kernel of pf,. Pi, is invariant 
under the transformation Tab = Ta ® Tb, where 

= ( J e^2*!/3 ) > (25) 

/I 0\ 

_ I cos(27r/3) -sin(27r/3) | 

^~ sin(27r/3) cos (27r/3) ' 

V 1 / 

Note that Tb is a real matrix. Later on we will need its 
eigenstates with real coefficients; they are |0) ± |3). Note 
also that „ = 1. 



B. Construction of nd— EW's 

We use now the methods developed in Section V to 
obtain a nd-OEW starting from pb- That is, we define 
Wb = Pi+ Pi , where Pi is the projector onto K{pb) = 
K[pI). Our procedure consists of first subtracting the 
identity to obtain Wi = Wb — eil. Then, we subtract 
P2 + Qi"; -F3 + QI"; stc. In the n-th step we will have 

Wn^Wr,-l~en{Pn+Ql), (26) 

where P„ (Qn) is the projector orthogonal to the space 
spanned by Pw^^i {Pw"^ We will use the symmetries 
of Pb to better understand the structure of Wn ■ 

(a) Wn = . We can prove this by induction. First, 

it is clear that Wi = Wi . Let us now assume that 
Wn^i = W^_i. Then we show that Wn = W^ . 
For that, we just have to show that the subspace 
spanned by Pw„^i is the same that the one spanned 
by Py^T , so that Qn = Pn- But this is clear since 

(b) TABWnT\g = Wn- Wc provc this by induction. 

First, for Wi ^ Pi -I- Pf - eil we have that 
TabWiTIs = TabPiT\b + TabPIt\j, ~ eil = 
Wu since TabPIt\j, = [TabPiT\bY (given the 
fact that Tb is real) and Pi is invariant under Tab - 
Then, let us assume that TABWn~iT\g — Wn-i- 
In order to show that TABWnT^^ — Wn we just 
have to show that Pn is invariant under Tab, or, 
equivalently, that the subspace spanned by Pw„-i 
is invariant under Tab- But this follows immedi- 
ately from the fact that TABWn-iT^^ — Wn-i- 
□ 

Starting the property (a) it follows that the vectors 
|e, /) e Pw„ will have |/) real (unless we have degen- 
eracies). This can be seen by noticing that those vec- 
tors minimize {e, f\Wn\e, f); defining We = (e|W„|e), we 
have that W^F = We — W\ is symmetric, and therefore 
the eigenstate corresponding to its minimum eigenvalue 
can be chosen to be real. On the other hand, starting 
from the property (b) it follows that if |e, /) S Pw„ then 
Tj^gje, /), T^^je, /) G Pwr,- According to that, we will 
typically have two kinds of product vectors in Pvi/„ : 

(1) |e, /) is an eigenstate of t\q with |/) real: There are 

only 4 possible product vectors which fulfill these 
conditions: {|0), |1)} ® {|0) |3), |0) - |3)}. 

(2) |e, /) is not an eigenstate of t\q-- Then, we will also 

have: T\^\eJ) and (T\sf\eJ) G Pw - 

We have carried out this procedure for pb and found 
nd-OEW for each 6. Wc find that for the optimal EW we 
have two vectors of the kind (1) and six of the kind (2). 
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In total we find eight product vectors in Pw , which span 
the whole Hilbert space and therefore the corresponding 
EW are optimal (see Corollary 2). This means that any 
operator of the form (||) with \ek, fk) G Pw the product 
vectors we have found, and > will be a full range 
separable density operator that lies on the boundary be- 
tween separable and PPTES. Up to our knowledge, this 
constitutes the first example of those operators . We 
have also created the PM corresponding to the nd-OEW, 
which are the first examples of non-decomposable PM 
with minimal "qubit" domain, or - equivalently - mini- 
mal hermitian conjugate codomain. 

In Fig. 1 we show for which h' pb' is still detected by 
the nd-OEW created out of ph. We find that for a given 
b, the optimal witness that we create detects all p^^ for 
b < b' . Thus, in the figure we plot b' as a function of 
6. As explained above, the corresponding positive map 
detects more than the witness itself. In the figure one 
can also see how much is detected by the positive map. 



xio 



*| POSIVITE MAP 



Witness 
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FIG. 2. Maximum A such that pt + Al is still detected by 
the witness and the positive map created starting from pi,. 
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FIG. 1. Values of b' for which if b < b' , is detected by 
the witness and the positive map created starting from pt- 

Obviously, the witnesses that we create do not only 
detect the density operators pb- For instance one can 
check how much one can add the identity to certain pb 
but still keeping the state entangled. That is for which 
A, Pb + Al is still detected by the witness. This is shown 
in the following figure. 



Finally, let us note that we have observed using numer- 
ical calculations that if one starts with a random projec- 
tor, P of rank 3, and optimizes the decomposable opera- 
tor = P -h in the same way as the one described 
here, then one will end up with a nd-OEW W, where p^^ 
is complete. This means that in order to create nd-OEW 
one does not need to know in practice an edge PPTES. 
In another words, optimization itself is a way to reach 
nondecomposableness. 



C. Analytical procedure 

In this subsection we will present an analytical way to 
create nd-EW's. Furthermore we will given an example 
of such a witness, which detects pb for all 6 G (0, 1). From 
Fig. 1 we see that the witness which detects most is the 
one we created out of pb, where b is very close to 1. We 
will work with the original pb (p^). 

We consider two hermitian operators A and B, with 
A positive on product vectors, i.e., {e, f\A\e, f) > 0, 
whereas B does not have to. As before we denote by 
Pa (Pb) the (not necessarily complete) set of prod- 
uct vectors on which A (B) vanishes. We require that 
for all |e,/) e Pa, {eJ\B\eJ) > 0. Then we define 
W{x) = ^{A + xB) for any real x. So we have the fol- 
lowing 

Lemma 9: If lim2;_,o(p)M/(a;) < then p is entangled. 

Proof: We prove that lim2,_o(e, /|IF(x)|e, /) > 0. 
This implies that if p is separable, then \\m.x^o{p)w(x ^ 
0. Let us therefore distinguish two cases: (i) if |e, /) € Pa 
then we have that hma;^o(e, /|II^(x)|e, /) = (e, f\B\e, /), 
which is, per assumption, positive, (ii) |e, /) ^ Pa then 
we have \\m^^Q{e, f\W {x)\e, f) = lim^^^of + where 
a = {eJ\A\eJ) > and 6 = {eJ\B\eJ). Thus this 
limit tends to infinity, which proves the statement. □ 

Note that W{x) is not an EW since it is not necessarily 
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positive on product vectors. However, one can make it 
positive by adding the identity operator to convert it into 
an EW. 

Corollary 6: Given any xq > 0, then W{xo) = 
^(A + xoB) + X^J , with A:,„ = - miniej) (e, f\^iA + 
XoB)\eJ) is an EW. 

Let us now iUustrate how we can use Lemma 9 to detect 
all the states pb- We define 
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(28) 



One can easily show that A — |'0)(V'l + {\'t'){4>\)'^^ : 
where IV') = |01) - |12) and |0) = |02) - |11). Thus 
this operator is positive on product vectors, since it is 
decomposable. Let us now use unnormalized states in 
order to present the set of product vectors on which A 
vanishes, i.e. Pa- Pa = Pai U Pas: where Pa^ — {(|0) + 
a|l)) (g) (a;|0) + 2/|3))Va, x, y} and Pa, = {(|0) + e**|l)) ® 
[x\0)+y\3) + z(|l) + e^**|2)]V$, x, y, z}. The operator B 
has to be positive on those product vectors, i.e., V|e, /) G 
Pa, (e, f\B\e, f) > 0. In order to show that this is indeed 
like that let us distinguish the two cases: |e, /) G Pai and 
|e, /) e Pa2- In the first case we have that 




(e|i?|e) = 



and so {e, f\B\eJ) 
\e,f) ePA2 then 



(elBle) 



-2a 





1 


-lap 


l + lap 










1 - 






-2a 





-2a* 


1 












= \x + y\^ - 




X 


-y\' 



(29) 



> 0. If 



/ 


2 


-2e** 





' \ 




-2e-** 


2 
















2 


-2e'* 


v 








-2e-'* 


2 / 



(30) 



which is a positive operator and so (e,/|P|e, /) > 0. 
So those two operators A and B fulfill all the required 
properties. Furthermore one can show that (p&)a = 
and {pb)B < for all < 6 < 1. Thus we have that 



\imx^o{W{x)pb) < for all < 5 < 1, where we defined 
W{x) = ^{A + xB). 

As mentioned above we can use now W{x) in order to 
create other PPTES just by adding product vectors on 
which W{x) vanishes. To find the product vectors we 
can add, all we need to do is to determine the intersec- 
tion between Pa and Pb- Since Pb = {(|0) + e'"*!!)) (g) 



[am 



-401 



\1) +b{\2) + e-"i'\3))]y(t),a,b} we have that 



S = PaCiPe 



Pa, HPb = {(|0) + e*^|l)) » (|0) 



e-'^\l} + 6-^2-^12) + e~*30|3^)V(/)}. Note that S spans 
a 5 dimensional subspace and that the orthogonal sub- 
space is spanned by the vectors {— 102) + |13),— 101) + 
|12),-|00) + |11)}. 



VIII. CONCLUSIONS 

Entanglement witnesses allow us to study the separa- 
bility properties of density operators. We have defined 
OEW, which are those that detect entanglement in an 
optimal way. We have given necessary and sufhcient con- 
ditions for an EW to be optimal, and we have shown a 
way to construct them. We have also concentrated on 
nd-EW, which are those that detect PPTES. We have 
extended the definitions of optimality and the optimiza- 
tion procedure to those EW. It turns out that one can 
optimize nd-EW by subtracting decomposable operators. 
We have also given an explicit method to construct nd- 
EW starting from "edge" PPTES. We have also men- 
tioned that this method works by starting out from ran- 
dom operators. We have extended our techniques to PM, 
and therefore given a method to systematically construct 
non-decomposable positive maps. We have illustrated 
our methods with a family of "edge" PPTES acting on 

(gff*. The corresponding PM constitute the first ex- 
amples of PM with minimal "qubit" domain, or - equiva- 
lently - minimal hermitian conjugate codomain. We have 
also constructed the first examples of separable states of 
full range that lie on the boundary between separable 
and PPTES. These states can be used for experimental 
realization of PPTES |3|. 

In this paper we have also introduced the "edge" 
PPTES, which violate the range criterion of separability. 
As shown in Appendix B, the "edge" PPTES allow us to 
construct a canonical form of PPTES in Hilbert spaces 
of arbitrary dimensions. They also allow us to give a 
novel sufhcient condition for non-separability which ap- 
plies to operators with positive partial transpose. It is 
based on the fact that among all PM (or EW) only the 
subset {Aedge} of those PM that detect edge PPTES are 
needed to study the separability of PPTES. This opens 
many interesting questions. Is it possible that in the set 
{Aedge} there is some map that is globally finer than the 
transposition? In another words, is there a map detect- 
ing the entanglement of all the states with non-positive 
partial transpose? What is the minimal subset of {A^dge} 
providing such condition? Is it finite? 
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Finally, let us consider the implications of the our 
results for the very interesting problem of locality of 
PPTES. There is a conjecture that those states can 
be local in the sense that they admit a local hidden vari- 
able (LHV) model for any set of possible local measure- 
ments. The problem is not trivial given the fact that it 
may be important to take into account the role of sequen- 
tial measurements and the possible existence of many 
copies. Quite recently it has been shown that PPTES 
satisfy Bell-type of inequalities introduced by Mermin 
pO| . It is not difficult to convince oneself that the set of 
states admitting LHV model for any fixed type of mea- 
surements is a convex set. Furthermore, extending the 
reasoning from it is easy to see that the set of sepa- 
rable states admits LHV models for any possible set of 
measurements. Hence, taking into account the results 
of this paper it follows that in order to prove, or to dis- 
prove locality of PPTES it is enough to study only "edge" 
PPTES. 

Note that the "edge" states have typically very small 
rank (the minimal rank is four in 3 (8) 3 systems, see Ref. 

There have been no examples of LHV models for 
states of low rank, so far. Thus, perhaps completely 
new techniques will be needed to study this problem. In 
this case the most symmetric PPTES provided recently 
[Q(c) seem to be the best suitable for the first test. 



In all this appendix we will use that given an EW, 
VF, and an operator P > we say that P cannot be 
subtracted from W if for all A > 0, — XP does not 
fulfill (I). In other words, there exist |e(A)) e Ha and 
|/(A)) e Hb such that 



(e(A),/(A)|(I^-AP)|e(A),/(A)) <0. 



(Al) 



Note that (e(A), /(A)|P|e(A), /(A)) must be strictly pos- 
itive, so that (Al) can be expressed as 



(e(A),/(A)|W-|e(A),/(A)) 
A^o (e(A),/(A)|F|e(A),/(A)) 



0. 



(A2) 



In the first subsection we will introduce some defini- 
tions and notation. In the second one we give a method 
to determine the set of product vectors Pwi on which 
W vanishes. In the third subsection we find a necessary 
and sufficient condition under which an operator cannot 
be subtracted from an EW. We will see that there must 
exist a vector |eo, /o) G Pw, some other vectors |ei) and 
l/i), and certain phases 4>ej and 9 such that some quan- 
tity is zero. In the next subsection we will see that the 
problem can be reduced to finding only the vectors |eo,i) 
and |/o,i). Finally, we will show that if dim(iJ/i) = 2 we 
just have to find |eo) and |/o), which is very simple. 
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APPENDIX A: OPTIMALITY OF EW 

In this Appendix we study necessary and sufficient con- 
ditions for an EW to be optimal. According to Theorem 
1 of Section III we have that an EW, VF, is optimal iff no 
positive operator can be subtracted from W while keep- 
ing property (I). This condition can be reexpressed in 
terms of the infimum of some scalar products in Lemma 
4. This infimum is, in general, difficult to calculate (at 
least analytically). In this Section we will give a different 
method to determine whether an EW is optimal or not. 
This method will turn out to be very simple for the case 
in which dim(if^) ~ 2. The idea is to find the condi- 
tions such that a given operator P > can/cannot be 
subtracted from an EW. This will give us automatically 
a criterion to determine when W is optimal. 



In order to prove the results of this appendix in a com- 
pact and readable form we have made an extensive num- 
bers of definitions. 

We will always denote by |eo,/o) a product vector in 
Pwi and by |ei) S Ha and |/i) e Hb two vectors or- 
thogonal to |eo) and |/o), respectively. We will use the 
following notation: 

Wlj = {e,J,\W\ekJi), {i,j,k,l = 0,1). (A3) 
and wc will write 



(A4a) 

'0,0 - 1^^0,0 1^ ■ (A4b) 
We will also define the following operators: 

(A5a) 
(A5b) 



^ {e.,\W\e,), 



The following vectors will be used in the context of Eq. 
1^ 



W)) 



1 



:(|eo)+ecos(^)e'^Mei)), (A6a) 



v/l + |cos(0)e|2 

1/(6)) = , (|/o)+esin(g)e-^/|/i)), (A6b) 

v/1 + |sm(6l)e|2 
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where e is a real number, and (pej € [0, tt) and 
G [0, 7r/2] are certain constants. Given a prod- 
uct vector |e(e),/(e)) and an operator, W, we will ex- 
pand (e(e), /(e)|M^|e(e), /(e)) by collecting terms with 
the same powers in e; that is, except for a normaliza- 
tion constant. 



{e{e)Jie)\W\eie)Jie)) (x ^ e^MW), (A7) 



where 
Ao{W) 



Ai{W) = 2Re 



A2{W) = cos^{9)W^;^ + sm^{e)Wl^;^ 
+2sm{6) cos{e) 

xRe [e-'('^<'-"^^)VF°;o +e*("^=+"^^%o;o 
A3{W) = 2sin(6') cos(6l) 
xRc 



(A8a) 
(A8b) 
(A8c) 



(A8d) 
(A8e) 



cos(6')e**^W^i;o + sin(6')e*"*=lFo 1^ 

On the other hand, we will define 

1*04) = sin(0)e^*^|eo,/i) + cos(0)e'^=|ei,/o). (A9) 

Finally, the following quantity will play an important 
role in determining whether there exist vectors and pa- 
rameters for which (A2): 



x{w) ^ <b°<i - (Ko I + Ko I) 



O.li 



rl4l\2 



(AlO) 



2. Determining Pw 

As stated in Lemma 3, not every positive operator P 
can be subtracted from an EW, W; it must vanish on 
Pw- Thus, in order to choose P one has to know the set 
P\Y- In this subsection we give a method to determine 
it. 

We start by characterizing the vectors in Pw- 
Lemma Al: Given an operator W satisfying (1), then 
|eo,/o) e Pw iff 



(eo|W^|eo)|/o) -0, 
(/o|l^|/o)|eo) =0, 



(Alia) 
(Allb) 



Proof: (If) We just apply (/o| to Eq. (|Alla| ). (Only 
if) Since W fulfills (1) then We,, = (en I I en ) mus t be 
positive. Thus, (/nlW'eol/o) ~ inr phes Eq. ( |Alla| ). In 
the same way we obtain Eq. ( [Alibi ). □ 

In practice, for a given W the set Pw can be found 
as follows. Due to the fact that W is an EW we have 
that for any |e) G Ha, We = (e|P^|e) must be a positive 



operator (i.e. {f\We\f) > for all |/) G Hb)- Thus, the 
determinant det(lVe) > 0. According to Lemma Al, this 
determinant is zero iff there exists some |/n) G Hb such 
that (/n|W^eol/o) = 0, i.e., if |eo,/n) G Pw- That is, the 
determinant as a function of |e) has a minimum (which 
is zero) at |en). We can use this fact to find |en). Then, 



we can easily obtain |/n) via Eq. (Alia). We can expand 
an unnormalized state |e) in an orthonormal basis {|fc)} 
as 



dim(_ffA) 

|e) = •='=1'^)' 

k=l 



(A12) 



and impose that the corresponding determinant is zero. 
This gives us a polynomial equation for the coefficients 
Cfe, i.e. 



P(cfc,4) = 0. 



(A13) 



We also impose that, given the fact that the determinant 
is a minimum, 

7^Pick,cl) = ^P{c,,4) = 0, (A14) 



dck 



del 



which also give a set of polynomial equations. These 
equations can be solved using the method mentioned in 
Ref. M. 



3. Necessary and sufficient conditions for 
subtracting an operator 

In this subsection we give a necessary and sufficient 
condition for an operator P to be subtractable from an 
EW. We start out by giving s ome properties of the coef- 
ficients A{W) defined above ( |A8|) . 

Lemma A2: Given W satisfying (1) and |en,/o) G 
Pw, then for all |ei) G Ha and |/i) G Hb we have 

(i) AoiW) ^ Ai{W) ^ 0. 

(ii) A2{W) > 0. 

(iii) If A2{W) = then ^3(1^) = 0. 

Proof: (i) It is a direct consequence from Lemma Al. 
In order to prove (ii-iii) we use the fact that W sat- 
isfies (I). We define |e(e)) and |/(e)) as in (^). We 
impose that(e(e), /(e)|VF|e(e), /(e)) > 0. Using the 
expansion (|A7|) and taking into account (i), we have 
A{e) = A2{W) + eA3{W) + €^Ai{W) > for aU e. This 
automatically implies (ii), since otherwise for sufficiently 
small e we would have A{e) < 0. It also implies (iii), 
since if AsiW) < (AsiW) > 0) then for sufficiently 
small e > (e < 0) we would have A(e) < 0. □ 

Now, we are at the position of giving a necessary and 
sufficient condition under which an operator cannot be 
subtracted from an EW: 



13 



Lemma A3: Given P fulfilling PPw = 0, it cannot be 

subtracted from W iff there exists |eo,/o) G Pw, |ei) -L 
|eo), l/i) -L l/o), 0e/, and 9 such that A2{W) = but 
A2{P) ^ 0. 

Proof: (If) We define |e(A)) and |/(A)) as in @. Us- 
ing Lemma A2(i) we have Ao{W) = Ao{P) = ^i(M^) = 
Ai{P) = 0. Using Lemma A2(iii) we have that A3IW) = 
0. Thus, we can write the limit (A2) as 



X^A^iW) 



x^o A2{P) + XAsiP) 



X^A^iP) 



(A15) 



which obviously tends to zero given that A2{P) 7^ 0. 
(Only if) There exist two normalized vectors |e(A)) and 
|/(A)) (continuous functions of A) fulfilling ([A^). Tak- 
ing the limit A — > in this expression wc have that 
(e(0),/(0)|W^|e(0),/(0)) = 0, and therefore |eo,/o) = 
|e(0), /(O)) £ Pw- This means that wc can always choose 
|e(A)) = |e[e(A)]) and |/(A)) = |/[e(A)]) given in (0), 
where |ei) _L |eo) and |/i) _L |/o) are two normahzed vec- 
tors, limA^oe(A) = 0, and (e(e), /(6)|P|e(6), /(e)) ^ 0. 
We use ( |A6D t o ex pand the numerator and denominator 
of (A2) as in (A_7). According to Lemma A2(i) we have 
that Ao{W) = AoiP) = Ai{W) = Ai{P) = 0. Thus, we 
must have 



A2{W) + eAsjW) + e^AijW) 
e"'5 A2{P) + eAsiP) + e^AiiP) 



lim 



0. 



(A16) 



This implies A2iW) = and A2{P) ^ 0. Note that if 
both ^2(1^) = A2{P) = then, according to Le mma 
A2(iii) we have that A3{W) = AsiP) = 0, so that ( |A15| ) 
would require A4{W)/A4{P) — 0. But this cannot be 
since A4(W^) = would imply that |e(e),/(e)) £ Pw, 
and therefore (e(e), /(e|P|e(e), /(e)) = 0. □ 

Finally, we show in the next lemma that condition 
A2{P) = is equivalent to having certain vector in the 
kernel of P. We will use the vector |^o.i) defined in (A£). 

Lemma A4: Given a positive operator P and a set 
of vectors |eo,/o) € K{P), |ei) _L |eo) |/i) -L |/o), and 
parameters (/)ej, and 9 then A2{P) = iff |*o,i) ^ ^i^)- 

Proof: Since P > and |eo, /o) G K{P) we have P°;° = 
0. Then, we can write A2{P) = (^'o,i|-P|^'o,i), with l^f) is 
defined in (A9), from which it is obvious that A2{P) = 
iff l^-oa) e K{P). □ 



4. Necessary and sufficient conditions for ^2(14^) — 

The previous lemmas tell us that we cannot subtract a 
given operator P provided we can find some vectors and 
parameters such that j42(M^) = 0. The task of finding 
these vectors is difficult, in general. Here we will give a 
way to check whether these vectors exist. As before, we 
will denote by |eo,/o) a vector in Pw, and by |ei) and 
1/1) two vectors orthogonal to the first two. The quantity 
X{W) defined in (AlO) will play an important role in 
determining whether there exist vectors and parameters 



for which A2{W) = 0. In this subsection, we will always 
have to choose the phases (pej that minimize A2{W). 
That is 



-i(0c-<^/-0o) _ _l gi(0e+0/ + 0l) 



-1. 



(A17) 



We will denote A2{W) the value of A2{W) for this par- 
ticular choice of phases. We have 



A2{W) = cos2(0)M^io° + sin2(0)M^o";i 



(A18) 



-2 sin(6l) cos(6') ^W^;°W°^^ - X{W), 



where we have used ( AlO[ ). 

Let us start showing that X{W) is positive. We will 
use this property later on to reexpress the condition 
^2(W^) = in terms of one that is simpler to check. 

Lemma A5: X{W) > 0. 

Proof: This follows from the fact that A2{W) > 
for all values of (t>e.f- In particular, A2(H^) > 0, which 



according to ([Al|j imphes X{W) > 0. □ 

The next lemma shows that we just have to check 
whether X{W) = if we want to see if there exist pa- 
rameters (j)ej and 9 such that A2(PF) — 0. This first 
condition is therefore much more useful than the last one. 

Lemma A6: X{W) = iff there exist and 9° 
such that ^2(1^) = 0. 

Proof: (If) Given the phase 9 = 9^ we have that = 
A2{W) > A2iW). Thus, A2{W) = 0. According to 
( [Al^ ) we can have two cases: (a) 9o 7^ 0, 7r/2. In that 
case it is obvious that X{W) = 0. (b) 9o = 0, n/2. In the 
first (second) case we must have W^'q = {Wq^ = 0). 



But this implies that Wr 



1,1 
0,0 



^o'l — since otherwise 



we could always find some other value of 9 such that 
A2{W) < 0. Then, X{W) = 0. (Only if) We choose (f>ej 
as in ( |A17| ). For this value, according to ( AI8 ) we have 



A2iW) 



cos{9)JW^;o - sm{9)JW^ 



/0,1 

0,1 



(A19) 



which can always be zero for some particular value of 9. 
□ 

Note that according to the proof of Lemma A6, if 
VF^'p" = then A2{W) = only for 9^0. But in that 
case one can easily check that the vector |e(A),/(A)) G 
Pw [see (|A6|)] which cannot be. Similarly, we conclude 
that W^i '0° 7^ if we want A2{W) = 0. Thus, fro m now 

one we will assume that both Wiq and W^q are not 
zero. 



5. Optimality test 

Thus, we can now state the steps to check whether an 
EW, W, can be optimized or not. (1) For each |eo, fo) G 
Pw we must check whether there exist |ei) _L |eo) and 
1/1) _L 1/0) such that X{W) = 0. Let us denote by \e^o\) 
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and I/qI) the set of vectors fulfilling that. (2) For each 
of these vectors, we have to find the corresponding val- 
ues of </)g*y by using ( A17 ) and of 0^*^ by imposing that 

A2iW) = in ( [ATg| ). (3) Construct |*(*)) according to 
([A9|). (4) See whether the space spanned by Pw and 
is equal to Ha®Hb- If it is, then W is optimal. 
If it is not, we can always find some \'ip) orthogonal to 
that subspace that can be subtracted from W . 



6. Necessary and sufficient conditions for X{W) — 

The hard part of the procedure outlined before to see 
whether and EW is optimal is the step (1), namely to 
find |ei) and such that X{W) — 0. We start out by 
giving a necessary and sufficient condition for X{W) = 0. 

Lemma A7: Given |eo, /o) G Pw^ and |ei) _L |eo) and 
|/i>±|/o>, then X(iy)-0 iff 



i.e. they are negative real numbers. We obtain that 

1 

"o.o 



e-*'*=(/oKo-^<ol/o) <0, 

Wn 



(A24) 



so that we determine cjie- With these results, we can 
prove the following necessary and sufficient condition for 
X{W) = when dim{HA) = 2. 

Lemma A8: If dim(iJyi) = 2, given |eo,/o) G Pw, 
then there exists |ei,/i) such that X{W) = iff 



(/ol 



1 



1 



0,0 



0,0 



l/o 



(/ol 



^o,i: 



-w'ail/o) 



1^0,0 



Proof: (If) We define 
1 



l/i 



-(e- 



-'1,0 ■ 



<i)l/o 



(A25) 



(A26) 



W, 



\ <o" 



(e-**=<o + e^^=<i)|/o). 



«^o,o|ei) 



where (j>. 
Proof. 



W, 



0,1 



\ <o" 



o,i)|eo), 



(A20a) ^ ^ _ 

where 0e is determined by the condition (A24). Using 
this expression to calculate X{W) one finds that indeed 
X{W) = . (Only if) Using Lemma A7 we can w rite 
(A20b)as in so that the phases 0e,/ ensure that ( |A23| ) is 

fulfilled. Substituting |/i) in the equation X(W) = one 



J are given in (A17). 
'(If) We multiply by 



Eq. ( |A20a| ) and take 
the square of the absolute value of the result. We obtain 



i,i| 
0,0 I 



\W 



i,'ol 



(A21) 



finds ( |A25| ). □ 

In summary, for a given |eo, /o) G Pw, in order to find 
whether there exist |ei, /i) su ch th at X(IU) = we just 
have to check the condition ( A25| ). If it is fu lfilled , we 
can e asily find |/i) and the phases (pe / using ( A23 ) and 



Using Lemma A5 we conclude that X{W) = 0. (Only if) 
Since X{W) = and according to Lemma A5 X{W) > 0, 
then X{W) must be a minimum with respect to |ei) and 
l/i). Taking the derivatives of X{W) with respect to 
these tw o ve ctors and imposing that they vanish, one 
obtains ( |A20| ). □ 

Equations (A20) are particularly useful if the dimen- 
sion of one of the Hilbert spaces is 2. Without loss of 
generality, let us assume that dim(iJyi) = 2. In that case 
we can choose |ei) as the one that is orthogonal to |eo) 
(with an arbitrary choice of the global phase) . The de- 
termination of (pe can be done as follows. Using (A2C) 
we write 



l/l) 



1 



w, 



ie-^'^^wlo + e^'^^wl,)\fo) 



0,0 



where 1 / wg q denotes the pseudo-inverse 
use this expression to impose 



(A22) 



We can 



(A23) 



APPENDIX B: CANONICAL FORM OF PPTES 

The concept of "edge" PPTES seems to play a very 
special role in the characterization of PPTES. In partic- 
ular, in view of the criterion given in Section VD, which 
is based on the fact that any density operator p can be 
decomposed into a separable part and an "edge" PPTES 
(|l8| ) . Among all the possible decompositions there might 
be one for which the trace of the separable part is max- 
imal. When it exists, such a decomposition was termed 
positive partial transpose best separable approximation 
(PPT BSA) to pJpL It extended the idea of BSA intro- 
duced in Refs. [ p3| , ^2[ to the case of PPTES, which were 
based on the method of diminishing the range of p by 
subtracting product vectors from its range, while keeping 
the remainder and, at the same time, its partial trans- 
pose, positive |23y22| , |20| , pT|| . In this Appendix we formal- 
ize the results regarding the existence and properties of 
the PPT BSA. In particular, the proofs presented in the 
quoted papers were restricted to the case in which there 
exist a finite, or at most, countable number of projectors 
on product vectors that can be subtracted from p. We 
will extend them below to continuous families of product 
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vectors. The Appendix is written in a self-contained way, 
and can be read independently of the body of the paper. 

We denote by Tp the set of projectors on product 
vectors {\ea, fa){ea, fa\} such that \eaj fa) G R{p) a-iid 
|e„,/*) e In Ref. (2l) we showed that if Tp is 

finite then there exist an optimal decomposition (PPT 
BSA) p = {l-p)psep +pS where 5 is an "edge" PPTES, 
and p is minimal. Note that PPT BSA involves the state 
6 which violates the range criterion in a rather special 
way, i.e. with the additional requirement that Tp is a 
finite set. It can happen that there is an uncountable 
family of product vectors depending on continuous pa- 
rameter that can be used for subtracting projectors. In 
the following we will show that in such case the above 
result is valid. 

In order to consider the case of continuous families of 
product vectors we first prove the following: 

Lemma Bl: Let p will be a PPTES defined on a 
Hilbert space Ti. , dimTi < oo. Then the set of product 
vectors Fp is compact. 

Proof : Obviously Tp is a bounded set in finite- 
dimensional space, so it is enough to show that it is 

e 



closed. Consider any sequence |(?„,/in 



R{p), \gn,K) € The limit vector must: (i) 

respect the condition of orthogonality to K{p) [i. e. 
they must belong to R{p)], (ii) belong to the sphere 
(i. e. set of all vectors with — 1), (iii) fi- 

nally, it must be a product state, because if it was entan- 
gled then its distance from the compact set of product 
pure states [|3| defined as min|e j) II If/*) — |e, /) || would be 
nonzero, which is obviously impossible. We conclude thus 
= \g, h) G R{p) for some \g), |/i), which implies (up to 
irrelevant phase factors) that \gn) — *■ \g) and |/i„) — )■ \h). 
We have (again up to irrelevant external phase factor) 
\gn,Ki) Is, /i*)- The latter must belong to R^p^'^) as 
any element of the corresponding sequence is orthogonal 
to K{p'^"). □ 

Let us now prove the following general lemma, which 
is a generalization of one theorem from Ref. : 

Lemma B2: Let the PPTES p be defined on a finite 
dimensional Hilbert space. Consider the set Sp consisting 
of the trivial zero operator plus all unnormalized states p 
(trp < 1) such that 6 = p — p\s positive and has positive 
partial transpose. Then, one can find p G Tip such that 
with tr(p) < 1 is optimal in the sense that: 

(i) The trace oi S = p — p is minimal with respect to all 

separable p's leading to positive partial transpose 

(ii) The state S = (5/tr(<5) is an "edge" PPTES. 

Proof : To prove the existence of p € Sp we just have 
to show that Sp is compact. This can be done by showing 
that Sp is a closed subset of another compact set, namely 
C — convlFpU }. The latter set C is compact as it is 
a convex hull of the compact set {FpU } in a finite 
dimensional space. 



Note first that Sp C C. Indeed, by virtue of 5 > any 
nonzero p cannot have any vector in its range not be- 
longing to R{p)- Analogously R{p^'^ ) c R{p^'^). Hence, 
according to the properties of the ranges of density oper- 
ators in general JlSf , p must be a convex combination of 
vectors from Fp, and as such it belongs to C. Let us show 
that Sp is closed. This follows immediately form the fact 
that Sp is a cross-section (performed over any projections 
P, Q) of the sets: S^ p = {p : fp^p) = tr{Pp~Pp) > 0} 
and S2 Q = {p : gQ,p{p) = tr(Q^«p - Q^«p) > 0}. Since 
the functions fp,p,gQ,p are continuous, all the sets par- 
ticipating in the cross section are closed. Now, the cross- 
section of closed sets is again a closed one. 

Consider now the statement (ii). Since (5,(5"^^ > 0, we 
always have 5 = /3Pr(5) -t- A and (5^« = P'Pr(si'b ) + A' 
with f3, f3' > 0, some positive operators A, A' (here, Px 
denotes a projector onto the subspace X C H). Then 
if, contrary to (ii), there were any |e, /) G R{S) such 
that |e, /*) e R{S'^^), then the new operator p* ^ p + 
7|e, /)(e, /I, 7 = min[/3,/3'] would fulfill that S* = p-p* is 
a PPTES, and would contradict optimality with respect 
to (i). □ 

Let us remark that if we give up the condition regard- 
ing positivity of S^'^ , then we obtain a modified state- 
ment (ii) where the state S has no product vectors in its 
range. This is nothing but the best separable approxi- 
mation (BSA) of Ref. |2^, extended here rigorously to 
the states p having uncountable set of product vectors in 
R{p). 

From the Lemma B2 we obtain the following charac- 
terization of PPTES, which can be regarded to be among 
the main results of this appendix, since it provides a 
canonical form of PPTES: 

Proposition : If the state p is PPTES, then it is a 
convex combination 



p = (1 -p)a 



Sep 



-pS (Bl) 

of some normalized separable psep and a normalized 
"edge" PPTES S. In the above decomposition the weight 
p is min ima l [i. e. there does not exist a decomposition 
of type ( Bl ) with a smaller p] . 

The above proposition means, in particular, that the 
edge PPTES are responsible for PPT type entanglement. 
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